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M 


otivated by recent reports of improved customer service facilitated by the use of transshipment to pool 
inventory, we develop a formal model that focuses on the role of transshipment in a system of n-retailers 


who stock a good. In addition, some stock may be kept at a warehouse. If a retailer faces a stockout, he makes an 
attempt to fill backorders by using a transshipment that may come either from the warehouse or from another 
retailer that has excess stock. We determine that the optimal transshipment policy is described by exactly one of 
five protocols, and that the choice among them can be made by evaluating a set of easily computable conditions. 
We also provide other structural results that help identify conditions under which a warehouse should be open: 
how this decision is influenced by the total number of retailers, and how this decision is related to the degree 


of correlation of demand between retailers. 


Key words: multilocation; transshipment policy; inventory pooling; demand correlation 
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Introduction 

In conventional models of stochastic inventory sys- 
tems a retailer has limited options when faced with an 
unanticipated stockout: backorder demand, if feasi- 
ble, or lose the sale. However, when operating within 
a network of facilities, opportunities are available for 
pooling inventory to improve performance. One such 
opportunity is to transship inventory on hand at one 
location to meet demand at another location that faces 
a stockout. 

The use of inventory pooling is quite prevalent and 
is increasingly supported by computerized inventory 
tracking systems. For example, a bookseller or cloth- 
ier will often identify another nearby retailer from 
the same chain that has the requested item in stock. 
Often, the item is put on hold and arrangements are 
made for customer pickup. An alternative to customer 
pickup is an expedited transshipment, as is the prac- 
tice for replacement parts in auto dealerships, from 
the service department of one dealer to that of another 
dealer that faces a stockout. 

Obtaining an expedited transshipment from 
another dealer is usually the first level of response 
to stockouts. If the service department is unable 
to locate a part from a nearby dealer, a warehouse 
or distribution center is used as a second level of 
response. Narus and Anderson (1996) describe sev- 
eral practical instances of the use of transshipments 
and multiple levels of responses in retail/ warehouse 
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networks. Depending on the protocol, the inventory 
tracking system may first search for the stocked out 
item at facilities at the retail level or the warehouse. 
In another protocol only the warehouse is a source 
for such transshipments, as is the case for the shoe 
company described by Retzlaff-Roberts et al. (1995). 
This is appropriate when organizational structures, 
like franchising arrangements, inhibit cooperation 
among retailers who are independent agents. 

It is the goal of this paper to study the perfor- 
mance of inventory systems with pooling in the fol- 
lowing one-warehouse-n-retailer system. Before the 
selling season, each facility places an order to bring its 
inventory on hand to its target level. After its random 
demand is observed, each retailer satisfies as much of 
its demand as possible from its inventory on hand. 
If economical, each retailer that has unfilled demand 
then makes an attempt to fill backorders from the des- 
ignated list of first-response facilities. If some demand 
remains unfilled an attempt is made to satisfy it from 
a second-response facility, if any are so designated. 

As we show in this paper, for our model the opti- 
mal policy must be one of the following five sys- 
tems. In the retailer-only system (RO) transshipments 
may only originate at the retail level. In the retailer- 
first system (RF) transshipments may originate from 
the warehouse only if all stock at the retail level has 
been exhausted. In the warehouse-first system (WF) 
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this protocol is reversed: Transshipments may origi- 
nate from the retail level only if inventory at the ware- 
house has been exhausted. In the warehouse-only sys- 
tem (WO), as the name implies, transshipments may 
not originate at the retail level. Finally, in the no- 
pooling system (NP) no transshipments are permit- 
ted so each retailer acts independently. Furthermore, 
we fully characterize the policy structure by precisely 
identifying the easily computable conditions under 
which each of the five systems are optimal. We also 
provide other structural results that help identify con- 
ditions under which a warehouse should be open: 
how this decision is influenced by the total number 
of retailers, and how this decision is related to the 
degree of correlation of demand between retailers. 

While related single-period models that incorpo- 
rate inventory transshipments to reallocate inventory 
have been extensively studied, most of the literature 
has focused on retailer-only systems. For example, 
in an early paper by Allen (1958), demand at each 
of n retailers is described by an independent nor- 
mally distributed random variable. At the beginning 
of the period a central decision maker has the oppor- 
tunity to redistribute initial inventories across retailers 
in anticipation of stochastic demand in an n-retailer 
network. At each retailer, realized demand in excess 
of the starting inventory is lost. Gross (1963) general- 
izes the model considered by Allen in two significant 
ways. In addition to the opportunity to redistribute 
initial inventories, there is now an opportunity to pur- 
chase additional units. Moreover, the demand dis- 
tributions need not be restricted to members of the 
normal family. While Gross assumes that all economic 
parameters are identical at each retailer, Karmarkar 
and Patel (1977) allow for differing values at each 
retailer. A related extension is due to Herer and Rashit 
(1999) who, for the two-retailer case, incorporate fixed 
charges for placing initial orders. 

In these papers all decisions are made before the 
start of the period—there is a single opportunity 
for outside procurement. This procurement is com- 
plicated by the reallocation decision. In contrast, 
Simpson (1959) proposes a different recourse strategy. 
In his model, a given quantity of inventory is to be 
allocated to n retailers at the start of the period, and 
hence the allocation cost is sunk. Then at the end 
of the period a recourse decision is made whence a 
retailer that faces a stockout has the option, for a fixed 
fee, of obtaining sufficient stock from an emergency 
supplier to meet backorders. Krishnan and Rao (1965), 
like Simpson, allow for a recourse decision at the end 
of the period. However, transshipments to a back- 
ordered retailer are made at constant unit cost from 
any retailer that has unused inventory. Tagaras (1989) 
extends this model to the case of nonidentical eco- 
nomic parameters but only considers the two-retailer 
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case. Recently Dong and Rudi (2004) have extended 
the model of Krishnan and Rao (1965) in two substan- 
tial ways. They allow the demand process at retailers 
to be drawn from an arbitrary multivariate normal 
distribution and provide some asymptotic properties 
of the optimal response for the case of identically 
distributed demand for a very specialized covariance 
structure. 

The papers reviewed above are a salient sample 
from the literature on RO systems. In contrast, sys- 
tems that hold reserve stock at a warehouse for satis- 
fying backordered demand at the retail level, as in our 
WE, WO, and RF systems, have received limited atten- 
tion. An important contribution is due to Hoadley 
and Heyman (1977). Their model formulation repre- 
sents the RO system as well as the two warehouse sys- 
tems, RF and WF. Under the assumption of identical 
cost parameters across and between retailers, Hoadley 
and Heyman prove the convexity of the cost func- 
tion. For their variant of the WF system, they also 
develop an easily computable condition for whether 
or not a warehouse should carry reserve stock. In 
addition, they demonstrate that for the RF system the 
equal fractile rule is optimal. This structural result for 
the RF system is analogous to the earlier result for 
the RO system due to Krishnan and Rao. In sharp 
contrast, there appears to be no structural analysis of 
the WO system. Subsequently, Rudi (2000) while con- 
sidering a dual sourcing problem used a modeling 
framework that subsumes the WO system as a special 
case. His results include demonstrating the convexity 
of the cost function, but in his case the demand pro- 
cess is multivariate. Recently, Herer et al. (2004) have 
developed a viable approach for computing optimal 
policies in these systems. 

We contribute to the literature on systems that pool 
inventory using transshipments by providing a com- 
prehensive analysis of the optimal structure of the 
transshipment policy. In §1 we present the general 
model. We show that an optimal system that incor- 
porates transshipment must be one of RO, WO, RF, 
and WF; otherwise, the NP system is optimal. Conse- 
quently, we limit attention to these four systems for 
inventory pooling. In §2, we synthesize and extend 
results on the RO system. These results play an impor- 
tant role in the subsequent analysis, in §3, of the 
three warehouse-based WO, RF, and WF systems. 
Finally, in §4, we close by identifying the scope of our 
findings. 


1. Model Formulation and 
Characterization of the 
Optimal Policy 


We consider the following one-warehouse-n-retailer 
inventory system. Let 1,1=0,...,, denote the index 
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for the warehouse and retailers 1 to n respectively. 
Before the selling period, the warehouse and the 
retailers place and receive orders from an outside sup- 
plier to bring their inventory up to S;, i=0,...,n. 
Each retailer then experiences random demand drawn 
from an arbitrary multivariate distribution. At the end 
of the period each retailer satisfies as much of its 
demand as possible from its inventory on hand. Sub- 
sequently, there is an option to satisfy any demand 
not met by a retailer’s on hand inventory, from on 
hand inventory at the warehouse, or from unused 
inventory at another retailer. Demand not satisfied by 
such transshipments is lost. The objective is to deter- 
mine the optimal stocking levels so as to maximize 
the expected profit (net of transshipment and penalty 
cost). In developing the model, we will use the fol- 
lowing notation. 


Cost Parameters 

c,: unit cost of understocking (lost sales) after all 
transshipments. 

C,/Cy: unit cost of overstocking at the retail level/ 
warehouse charged on ending inventory on hand 
after all transshipments. 

c,/Cy: unit transshipment cost from the retail level/ 
warehouse, including a late penalty cost if any. 

We assume that all parameters are positive. The 
parameters c, and c, with finite values are required 
for an “independent” newsvendor problem. To admit 
transshipment at the retail level, c, must be finite. To 
admit transshipment from the warehouse, c;, and c, 
must be finite. As we show later, the relative values 
of c,, cr, and c,, to those of c, and c, determine which 
transshipment protocol is optimal. In general, the unit 
cost of overstocking includes the physical storage cost 
and capital cost (including the unit purchase cost) 
while the unit cost of understocking includes loss 
of goodwill and unit contribution less unit purchase 
cost. In addition, as in Herer and Rashit (1999) the 
unit transshipment cost may include a late penalty 
cost that represents compensation to the customer for 
accepting a transshipment. 


Decision Variables 

Sp: stocking level at the warehouse at the beginning 
of the period. 

S;: stocking level at retailer 1, 1=1,...,n, at the 
beginning of the period. 

S: total system stock at the beginning of the period; 
S = Vino 5: 

S°/S;: optimal total system stock/stocking level at 
the warehouse or at retailer i, i= 1,...,n, at the 
beginning of the period for system *, e = NP, RO, 
WO, RF, WE. 


The Retail Level 
D: vector of random demand at each retailer; 
D= (Dj) 


Sa n° 
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4, o: mean and standard deviation vectors of D 
respectively; p = (14;);, © = (9;);- 

2: covariance matrix of D; %, = (0;;);,; where oj = 
p,,0;0; for iF]. 

F./f;; marginal cumulative distribution function 
(c.d.f.)/marginal probability density function (p.d.f.) 
of D,. 

F: joint cumulative distribution function of D. 

L;(S;): excess inventory at retailer 1 after demand is 
satisfied by inventory on hand. 

Lr(So,--+,5,): excess inventory at the retail level 
after all transshipments. 

Tr(So,.++,5,): total number of transshipments orig- 
inating at the retail level. 


The Warehouse 

Do(S;,.--,5,): random demand at the warehouse. 

H,,/h,(S,,...,5,): cumulative distribution func- 
tion/probability density function of Do, the demand 
at the warehouse when there are n retailers in the sys- 
tem; note that there is a mass at D) = 0. 

Lw(So,...,5,): excess inventory at the warehouse 
after transshipment from the warehouse to all stocked 
out retailers. 

Twr(So,--.,5,): total number of transshipments 
originating at the warehouse. 


The System 

fu: average aggregate demand in the system; 
b= Vij Mi- 

G,/g,: cumulative distribution function/proba- 
bility density function of random aggregate demand 
when there are n retailers in the system. 

BO(S,,.-.,5,): unsatisfied demand at the end of the 
period after all transshipments. 

L(S,...,5,): total excess inventory in the system 

after all transshipments. 
Note that BO, Tp, Twr, Lr, and Ly are random quan- 
tities that are linked to the random demand. It is also 
worth mentioning that Le, Lwr, BO, Tr, and Typ are 
the aggregate performance measures that will be used 
in computing the expected total cost of operating the 
system, which in turn is used to give the expected 
profit. 


1.1. Model Formulation and Structure 

Let E(-) and P(-) denote the expectation and the prob- 
ability of occurrence of “-” respectively. Our inventory 
model is a two-stage stochastic optimization prob- 
lem with recourse. In the first stage S,, the order- 
up-to level at each of n+1 locations is determined 
while in the second stage, after demand is real- 
ized the recourse decision of how to make trans- 
shipments to meet unfilled demand is made. As in 
the standard newsvendor model, we can write the 
expected profit function as expected riskless contri- 
bution less E[C(S), S,,...,5,)], the expected cost that 
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occurs as a result of the presence of uncertainty and 
includes an appropriate charge of lost contribution. 
Since the expected riskless contribution is indepen- 
dent of the decision variables, maximizing expected 
profit is equivalent to minimizing E(C). Let z; = S,—D,; 
denote the inventory level at retailer i after it has filled 
as much of its own demand as possible from inven- 
tory on hand; then z; denotes the leftover inventory 
on hand and z; denotes residual demand. Since the 
cost parameters are identical across retailers, when 
measuring the performance at the retail level, it suf- 
fices to use aggregate flows. Hence, the formulation 
of the mathematical program is as follows: 


min E(C(S, $1, ++, 8,)) 
=min | V(Sp,S, —D,,--., S, —D,) 0F(D), 
5S; YD 


where 
Vos Zip ewer Sa) 


= min C, Tp + Crlwr + Cy, (x ra — Ts | 
Tr, Twr ; 


+ Cy(Sp — Twr) + ¢, (x Zi — Iwr- Ts | / 
(P) 


subject to Tyr < 5p, 
Ip<) 2}, 
Twr + Ip eae and 


Twr> Tr = 0. 


To maintain focus on the strategic use of the ware- 
house for transshipment to meet unfilled demand, 
we explicitly exclude the option of “parking” leftover 
inventory at the warehouse at the retail level; and, 
conversely, “parking” leftover inventory at the retail 
level at the warehouse. It is easy to see that in a more 
general version of (P) that explicitly allowed “park- 
ing,” it would not be invoked in an optimal solution 
if the following two parking conditions hold: 


(PC1) 
(PC2) 


Cy —Cy SC; 
Cy — Cy < Cr. 
Hence, we shall assume that these two conditions 
hold. Consequently, it follows that 
THEOREM 1. If the condition 


(CO) 


Cr SC, 


does not hold, then in an optimal solution So is zero. 


Proor. The proof for this intuitive result is by con- 
tradiction. Suppose the optimal solution has positive 
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stock at the warehouse. There are two cases. If the 
demand realization is such that the optimal value of 
Twr is zero, reducing S, to zero saves cS). Alterna- 
tively, if Twp is positive, reducing it to zero saves the 
additional amount (cr; —Cp)Tyr. O 

To obtain additional insight into the optimal policy 
structure we now consider the second-stage problem 
in program (P). Using the fact that the second-stage 
problem in (P) is a transportation problem with two 
origins and one destination, it is easy to show that 
the optimal policy is quite appealing: Satisfy as much 
residual demand as possible, first depleting all stock 
from one level of the network and then using stock 
from the other level. To see this, first note that the 
objective function in the second stage formulation 
in (P) is equivalent to: 

min (C; — Cy — C,) Tp + (Cr — CH — Cy) Tyr- (1.1) 
Rs “WR 

Moreover, ruling out parking inventory in (P) leads 
to the following: 


T, < min} Y>2F, 32 - Twn} 


(1.2) 


(eee min{ Sy, ae Tat (1.3) 


To determine whether a facility will be used as a 
transshipment node and whether unsatisfied depend 
would be first satisfied by the warehouse or other 
retailers, consider the following three conditions that 
are derived by rearranging the coefficients in (1.1): 


(C1) 
(C2) 
(C3) 


C=) GC, 


Cr—CySc,, and 


Cr — Cy S Cy — Cy. 


Since the objective is to minimize the value of (1.1), 
it follows that if (C3) holds, using the warehouse as 
the first choice for satisfying residual demand would 
be the optimal response; otherwise, using the retailer 
first would be the optimal response. Moreover, from 
the linearity of (1.1), it follows that if a facility is used 
for transshipment, as much as possible of the resid- 
ual demand is satisfied from it. This rules out par- 
tial transshipment as an optimal response, allowing 
us to limit attention to the five systems, NP, RO, WO, 
WE, and RE. To choose between these five systems, 
observe that if (C1) does not hold, T, = 0, so that 
the RO and RF systems cannot be optimal. Similarly, 
if (C2) does not hold, Twr =0, so that using the WO 
and WEF systems cannot be optimal. If neither (C1) 
nor (C2) holds, using transshipments is not optimal. 
These observations are summarized as: 


LEMMA 1. In the optimal inventory system, the 
response to residual demand using transshipment satisfies 
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the following: 

(i) Partial fulfillment of residual demand (defined as 

above) using transshipment is not allowed. In particular, 
(a) if (C1) holds, then BOL, =0, and, 
(b) if (C2) holds, then BOLy =0. 

(ii) The response is hierarchical in that stock from one 
level of the network (warehouse or the retail level) will have 
to be depleted before stock from the other level is used. In 
particular, a transshipment from the warehouse is preferred 
to one from the retail level if and only if (C3) holds. 


Because (C0) implies (C2), it follows from Theo- 
rem 1 and Lemma 1 that the optimal system structure 
is given by one of the following five systems. 


THEOREM 2. In an optimal solution, the optimal sys- 
tem is one of the NP, RO, WF, WO, or RF systems. In 
particular: 

(1) if both (CO) and (C1) do not hold, then the optimal 
solution 1s given by the NP system; 

(2) if (C1) holds but (CO) does not, then the optimal 
solution is given by the RO system; 

(3) if (CO) holds but (C1) does not, then the optimal 
solution is given by either the WO or NP system; and, 

(4) if both (C1) and (CO) hold, then 

(a) if (C3) holds, then the optimal solution is given 
by either the WF or RO system, and, 

(b) if (C3) does not hold the optimal solution is given 
by either the RF or RO system. 


Theorem 2 provides a compelling result: By simply 
examining the cost parameters, one can identify the 
(only) benchmark system (where the warehouse is not 
utilized, i.e., the NP or RO system) and the (only) can- 
didate warehouse system (the WO, WE, or RF system) 
that need to be considered for the optimal inventory 
system. Another strategic implication of Theorem 2 
follows from reinterpreting (CO), (C1), and (C2) in 
terms of cp. Theorem 2 then reveals that the optimal 
transshipment policy depends crucially on the relative 
value of cp, as summarized in the following corollary: 


COROLLARY 1. AS Cp increases: 

(a) if (C3) holds and c, —c, < cr, the optimal policy 
evolves according to Figure 1(a); 

(b) if (C3) holds and c, —c, > cr, the optimal policy 
evolves according to Figure 1(b); and, 

(c) if (C3) is violated, then the optimal policy evolves 
according to Figure 1(c). 


Several interesting implications can be drawn 
from Figure 1. First, consistent with Theorem 1, a 
warehouse-based system is never optimal for c, < cr. 
More important, the WO system plays a limited role 
since it is admissible only if cr <c,—c, and only for 
intermediate values of c, between cr and c, — c,. For 
larger values of c,, the optimal system is one of the 
RO, RF, and WF systems. In view of Figure 1, in §3 we 
conduct a more thorough analysis of the effect of cp 
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Figure 1 Optimal System as c, Varies 
NP | NP | RO | RO or WF 
Cr—- CH Cr— Cp CT Cp 
(a) (C3) holds and c,—c), S$ cr 
NP | NP | NP or WO | RO or WF 
CT _ CH CT Cy _ Ch Cp 
(b) (C3) holds and c, -—c, > cr 
NP RO RO RO or RF 
Cr— Cp Cr— CH CT Cp 


(c) (C3) is violated 


on the optimal system. In particular, under suitable 
conditions, we find that for any given set of system 
parameters, there is a threshold such that for any c, 
larger than the threshold, the optimal system is either 
the RF or WF system depending on whether or not 
(C3) holds. 

While Theorem 2 and Corollary 1 are powerful, 
they do not provide a complete characterization of 
the optimal inventory system since in some cases 
(Cp > max{cr,c, — c,}) the candidate optimal inven- 
tory system is only reduced to two systems. However, 
because an important consequence of Theorem 2 and 
Corollary 1 is that in determining the optimal inven- 
tory system we can limit attention to one of the five 
generic systems, the objective function E(C) in prob- 
lem (P) can be simplified to 


E(C) = f (c,Lp + ¢,BO + ¢,Tp + CyLy + Cr Tyr} 9F(D) 


+ Cy E(Lw) + crE(Tyr)- (1.4) 
By adapting the standard result from the no pool- 
ing system that expected leftovers and backorders are 
convex functions of the stocking decisions, (1.4) can 
be easily manipulated to show that in the admissi- 
ble parameter space for each of the five systems, pro- 
gram (P) is a convex programming problem. 


THEOREM 3. 
(1) The expected cost function E(C) for the NP system 
iS convex. 
(2) If (C1) holds, then the expected cost function E(C) 
for the RO system is convex. 
(3) If (C2) holds, then the expected cost function E(C) 
for the WO system is convex. 
(4) If both (C1) and (C2) hold, then 
(a) if (C3) holds, then the expected cost function E(C) 
for the WF system is convex; 
(b) if (C3) does not hold, then the expected cost func- 
tion E(C) for the RF system is convex. 
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Table 1 K-K-T Conditions for the Five Systems 
NP RO 
(c, + ¢,)F(S;) —¢, — A; =0, /=1,...,2. 1.5) (C, +0, — C,)G,(S) + ¢,F,(S;) —¢, — A; =9, i=1,...,7. 1.8) 
A,$,=0, f=1,...,9. 1.6) ApS,=0, f= 1jpcwegll (1.9) 
A;,§;>0, #=0,...,A. (1.7) Aj,§;>0, f=1,...,0. (1.10) 
WO WF 
(Cy +0, — Cr)H, (So) + Cr — C, — Ag = 0. (1.11) (C, +C, — C,)G,(S) + (Cy + C, — C, — C7)H, (Sp) + C7 — C, — Ay = 0. (1.16) 
(Cy + Cp)F)(S;) + (Cy + 6, — Cr)dE (Ly)/AS; — c, — 2 =0, (Cy + Cy — C)G,(S) + (Cy + C — Cy — Cr)dE (Ly)/AS; + 6,F;(S;) — C, — A; =0, 
(27 cradle (112) j=1,...,n. (1.17) 
Ap Sp = 0 (1.13) Ap Sy = 0. (1.18) 
Apop=0). P= Misc (1.14) Apo; 0, fe 1 ycacg ll: (1.19) 
A;, 5; = 0, i=0Q,..., (1.15) A;, §; = 0, i=0,...,7. (1.20) 
RF 
(Cy +0, —C7)G,(S) +C; — C, — Ay = 0. (1.21) 
(Cy +0, —C7)G,(S) + (C, + Cr — Cy — C,)G, (x si) +¢,F,(S;)—¢,-—A,;=0, /=1,...,7. (1.22) 
i=1 
Ap Sy = 0. (1.23) 
A,S;=0, §=1,..., 0. (1.24) 
4;, 8; = 9, /=0,...,7. (1.25) 


It is noteworthy that the convexity result is based 
on condition (C2) and is thus slightly less restric- 
tive than Theorem 2 that uses condition (C0). Con- 
sequently, in the range cy — Cy < ¢, < Cr, while the 
relevant warehouse system has a convex objective 
function, using the warehouse can not be optimal. 

The most important significance of the convexity 
properties presented in Theorem 3 is that the Karush- 
Kuhn-Tucker (K-K-T) conditions, presented in Table 1, 
are necessary and sufficient to characterize an opti- 
mal policy. It is useful to observe that for systems NP 
and RO, because we assume that Ci O, it is never 
optimal to have zero stock at any retailer. In light of 
Theorems 2 and 3, we shall assume henceforth in our 
analyses that (C1), (C0), and (C3) apply when con- 
sidering the WO system; (C1) applies when consider- 
ing the RO system, (C1), (C0), and (C3) apply when 
considering the WF system; and (C1), (CO), and (C3) 
apply when considering the RF system. 

To proceed with the development of the three 
easily verifiable necessary and sufficient conditions 


presented in Table 2, which we name (NSyo), (NSwr), 
and (NSp-) respectively, under which it is optimal not 
to carry stock at the warehouse for each of the three 
warehouse-based generic systems, we need the fol- 
lowing characterization of dE(Ly)/dS; in the case of 
the WO or the WF system: 


LEMMA 2. In either the WO or WE system, dE(Ly)/ 
dS; = P({D; > S;} {Do < Sp})- 

Now consider the case when it is optimal to have 
Sy = 0 so that the warehouse is closed. When this 
is the case, it follows from Lemma 2 that dE(Ly)/ 
0S, =0. Moreover, when the warehouse is closed, the 
optimal response is that of the benchmark system 
since any inventory beyond that has no economic 
value because of the no-parking condition (PC2). 
Consequently, direct substitution into the appropri- 
ate K-K-T conditions followed by some algebraic 
manipulation results in the conditions in Table 2. An 
immediate consequence of Table 2 leads to a refine- 
ment of Corollary 1. Direct comparison of @wo with 
@yr reveals that they are equal when c, = c, — Cy, 


Table 2 Necessary and Sufficient Conditions Under Which the Warehouse Is Closed 
Owo 29, (NSyo) Owe = 9, (NSwe) Opp 20, (NSpr) 
where where where 
_ C(ONP npy _ __ &p — ©F _ F(ghO RO\ _ c,F,(S}") Cr _ lply + (Cp —Cr)(C,— Cy) F.(gR0 
Ovo =F (Sy, ...,8,° ) ray ey Owe =F (Sy, ....5,) Cy —C, +0, —Cr Gc=6.50 —t, RFC (Cp + Cy—Cr) (54). 
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leading to: 


COROLLARY 2. In the neighborhood of c, = Cc, — Cj, if it 
is optimal to use the WO system for c, < C,— Ch, then it is 
optimal to use the WF system for c, > Cc, — Ch. 


The significance of this result, as will be seen in 
Corollary 5, is that under appropriate regularity con- 
ditions it helps refine the conclusions of Corollary 1. 

We close this section by noting that our results 
demonstrate that (i) examining the economic parame- 
ters alone yields the optimal inventory system for the 
cases in parts (1) and (2) of Theorem 2, and (ii) exam- 
ining the economic parameters and the correspond- 
ing easily verifiable necessary and sufficient condition 
yields the optimal inventory system for each of the 
remaining cases. Moreover, determining whether or 
not a warehouse is open depends solely on the eval- 
uation of terms that depend on the solution to the 
NP or RO system. While the NP system corresponds 
to solving an independent newsvendor problem for 
each of the n retailers, the solution to the RO system 
is more complex and pivotal to our analysis. Hence, 
we now present a thorough analysis of this important 
system. 


2. The Analysis of the Retailer 
Only System 


In this section we focus on characterizing the stock- 
ing decisions S*°, the optimal stocking levels in the 
RO system. Since c, is positive, each retailer carries 
positive stock in an optimal solution. Therefore (1.8) 
reduces to 


(c, +c, ~¢,)G,(S*°) +¢,F(S°°)—cp=0, i=1,...,n. 

(2.1) 
From (2.1), as in Krishnan and Rao (1965), it follows 
that F,(S%°) = E,(SR°) =... = F,(S®°). Hence, the equal 
fractile rule is optimal. Without loss of generality, let 
F,(S%°) denote the value of the optimal fractile at 
each retailer. Although the computation for the opti- 
mal solution is reduced to a trivial search for one 
unknown, it is not possible to write this fractile inde- 
pendent of the demand distribution. This is because 
F,(SR°) and G,,(S®°) interact in (2.1). In contrast, the 
solution for the NP system is absolutely simple since 
it follows from (1.5) that each SN”, the optimal order 
quantity in the NP system, satisfies F(S‘") = c,/ 
(c, + c,), which only depends on each retailer’s 
marginal distribution. Consequently, the answer to 
the natural question: How does inventory pooling affect 
the optimal fractile relative to the NP system? requires 
significant effort. We proceed by first developing the 
following. 


THEOREM 4. 
(i) G,(S*°) = c,/(cy + ¢,) > F(Sf°) = c,/(c, + Cy) = 
F,(S)"). 
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(ii) Gic) < Col (Cy a Cn) ? F,(Si°) - Cy/ (Cy — Cn) vs 
F,(S)"). 

(iii) G,(S®°) > c,/(c, + cy) > F(SP°) < ¢,/(Cy +c) = 
F,(S}"). 


Proor. The theorem follows by substituting the 
condition in each case into (2.1). O 

In (i), each retailer’s service level and the system’s 
service level are equal to the newsvendor fractile. 
In (ii), each retailer’s service level is greater than the 
newsvendor fractile while the system’s service level 
is less than the newsvendor fractile. Finally in (iii), 
the inequalities in (11) are reversed. The key to opera- 
tionalizing Theorem 4 is to notice that if the marginal 
distributions and therefore the joint distribution are 
members of the normal family, a natural symmetry 
arises at the pivot point c, = c;,: 


COROLLARY 3. Assume that demand has a multivariate 
normal distribution. Then 

(i) If FS) = c,/(cy + Cy) = 1/2 (or c, =c;), then 
G,,(S®) = F, (SR) = 1/2, 

(ii) If R(SS") = c,/(c, + Cy) < 1/2 (or c, < c,), then 
G,(S®°) < (SN?) < B(SR°) < 1/2, and 

(ii) If Fy(SN*) = o/(cp + &y) > 1/2 (or c 
GG") > FS, Se (GS 2, 


» > Cy), then 


Hence, transshipping inventory maintains the same 
inventory investment if the newsvendor fractile is 
equal to half, increases the investment if the newsven- 
dor fractile is less than half, and reduces the invest- 
ment if it is greater than half. The precision provided 
by Corollary 3 leads to a full characterization of how 
the optimal response is influenced by variability: 


THEOREM 5. Assume that demand has a multivariate 
normal distribution with covariance matrix, %,. Then 
F,(S%°) is lexicographically increasing (decreasing) in x, 
while G,,(S®°) is lexicographically decreasing (increasing) 
Wt 2g f CC, (= C,). 


Theorem 5 provides the intuitive result that as 
variability increases the benefit of pooling decreases 
because each retailer orders more, and, at the same 
time, the probability that the system can meet 
demand falls. Analogously to Theorem 5, one would 
expect the same effect as the number of retailers, n, 
increases, with the benefit of pooling increasing with 
increasing n. While we are unable to provide a full 
characterization as a function of n, the analysis of the 
following two special cases shows that this insight has 
to be qualified. 


THEOREM 6. Assume that demand has a multivariate 
normal distribution. Suppose uw; = uw and o,; =o for 
all i, and oj = po*, p = 0, for all i # j. Then as 
n increases, F,(S*°) decreases (increases) while G,,(S®°) 


increases (decreases) if C, > Cy (Cy < Cj). 
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While it turns out that it is indeed the case that 
the benefit of pooling increases with increasing n 
when c, > c;, for the case when c, < c,, the oppo- 
site effect prevails, which is the consequence of the 
reflexive property of Theorem 4. As the next special 
case shows, this qualification to the benefit of pooling 
is not limited to the special case of identical retail- 
ers with a restrictive covariance structure; it can be 


strengthened for the case of independent retailers. 


THEOREM 7. Suppose demand at each retailer is inde- 
pendent and has an identical normal distribution. Then, 
(i) If c, < c,, then asn—o, 

(a) if c, = c,/2 so that c,/c, = 1/2, then F(S°) 
increases to 1/2 and G,(S®°) decreases to (cp — c,/2)/ 
(Cy, + Cp — C,); 

(b) if c, <c,/2 so that c,/c, < 1/2 and (cp — c,/2)/ 
(Cc, + Cp — C,) < 0, then F,(SK°) increases to C,/C, and 
G,,(S®°) decreases to 0. 

(ii) If Cy Cis then as n — oo, 

(a) if c, < 2c, so that (c, — c,)/c, < 1/2 and (c, — 
c,/2)/(Cy +¢, — C,) <1, then F,(SK°) decreases to 1/2 and 
G,,(S*°) increases to (c, — ¢,/2)/(c, +¢, — ¢); 

(b) if c, > 2c, so that (c, — c,)/c, > 1/2 and (c, — 
C,/2)/(Ch + Cy — c,) > 1, then F,(S*°) decreases to (c, — 
c,)/c, and G,,(S®°) increases to 1. 


While it follows from Theorem 6 that lim,,, ., F,(S®°) 
and lim, ,,. G,(S*°) exist when demand at each re- 
tailer has an identical normal distribution with p > 0, 
these limits can be identified exactly in the indepen- 
dent case so that they generate much more insight 
into the asymptotic behavior of the RO system, as 
well as the RF and WF systems as we will show 
later. For example, in (i)(a) and (ii)(a), S° converges 
to #4; aS n — co and hence, S®° converges to py, 
but G,(S*°) need not converge to 1/2. A possible 
explanation is the following: Let S, be an optimal 
stocking level at each (identical) retailer. The (com- 
mon) safety factor of each retailer and the system 
safety factor are defined as z, = (S, — w,)/o, and 
z= (nS, — np,)/./no, = Z,/n respectively. Note that 
although z, converges to 0 as n > oo, ./n diverges 
to co so that z = z,,/n need not converge to 0 and 
hence, G,,(S*°) need not converge to 1/2. It turns out 
that the rate of convergence of z, is proportional to 
/n so that G,,(S*°) converges to a finite positive num- 
ber, (c, — c,/2)/(c, +¢, — ¢,). On the other hand, (i)(b) 
and (ii)(b) are intuitive. In (i)(b), F,(S%°) converges to 
c,/C; (<1/2) instead, thus z, converges to a finite neg- 
ative number and hence, z = z,,/n diverges to —oo 
so that G,,(S*°) converges to 0. In other words, each 
retailer stocks less than its mean demand in the limit 
so that the probability that all system demand can 
be met is zero, ie., G,(S*°) = 0. Similarly, in (ii)(b) 
each retailer stocks more than its mean demand in the 
limit so that with probability one all system demand 


can be met, i.e., G,(S%°) = 1. In this case, we have 
a nice newsvendor interpretation since in the limit 
demand is always satisfied, either immediately or by 
delayed delivery. Hence, for each retailer, the cost of 
understocking is c, —c;, (the unit delayed delivery cost 
less the unit inventory holding cost at the source of 
delayed delivery) while the cost of overstocking is c,. 
Thus F,(S%°) approaches (c, — c,)/c, in the limit as 
nN —> oo. 

It is important to point out that the structural 
results on F,(S*°) presented in Corollary 3 and Theo- 
rems 5 and 6 are due to Dong and Rudi (2004). Apply- 
ing Theorem 4, we are able to extend these results by 
incorporating analogous results for G,,(S®°). We close 
this section by pointing out that one consequence 
of Theorems 6 and 7 is that the reflexive properties 
in n that are identified complicate significantly the 
characterization of when a warehouse is open in WF 
and RF systems, but the characterization of when the 
warehouse is open in WO systems is not affected since 
it uses properties of the NP system whose optimal 
inventory policy is independent of n. 


3. The Analysis of the WO, RE, and 
WE Systems 


In this section we will use properties of the bench- 
mark NP and RO systems to examine the strategic 
decision of whether to use a warehouse for satisfy- 
ing unfilled demand in the WO, RF, and NP systems. 
Since its resolution depends on the structure of ©, we 
first focus on developing its properties. To intuit the 
structure that we are seeking, consider the important 
case of n=1 for which the WO, WF, and RF systems 
are equivalent. Then ©, because of (1.5), simplifies to 


Cyl (Cy + €},) — (Cy — Cr) /(Cy + Cy — Cy) 
= [C,Cr — C,(Cy — CH) I/[(C, + Cn) (Cp + CH — Cr) I, 


which can be shown to be quasiconvex in c,. 
Since (C0) holds, the denominator is positive so that 
the sign of © depends on the sign of the numera- 
tor. Hence, © is positive when c, =0, and approaches 
zero aS Cc, approaches ov. If c, < cy, ® decreases to 
Cy —C}, (20) as c, increases. Conversely, if c;, > cy, © is 
quasiconvex and approaches c;, — c, < 0 from below 
as c, increases. Therefore, it follows that there is a 
threshold value of c, at which © becomes negative 
and remains negative for all larger values of c, if and 
only if c, > cy. This leads to the following: 


THEOREM 8. If n=1, SX° > 0 Sc, > cy and cp > 
Cyr / (Cy — Cy): 

This simple theorem has significant implications. 
First, notice that because there is only one retailer the 
motivation for holding stock at the warehouse can- 
not be based on an inventory-pooling argument, but 
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Table 3 Characterization of 0. as a Function of >,, , cp 
WO 
Condition for Cy >Cy/n (=) 
lM, 200 9 < OF 
Cy Quasiconvex 
(Ce or AN***) 
2. Cp <C, 
Co>C, + (AN) 
n Cp <C, 
Co >C, + (A) 
Necessary and (i) | 
sufficient (ii) Cp > Cy 
conditions for 
lim, 8 <0- 


N— oo 


Decreasing (A****) 


(i) 1.1.NP 
(ii) €,/2 <C,—Cy <C; 


(iii) c, > 
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RF WF 


Ch>Cy () Ch > Cy (<=) 


Quasiconvex (AN) 


t (AN) + (AN) 
1 (AN) Indeterminate 
J (1.1N) J (1.1N) 
* (LN) Indeterminate 


(i) I.I.N and one of 
the followings: 
(ii) ¢; <¢,/2 


2C,C; — C;(C; + Cy) 
ill) (C+ C7)/2 < 04/2 <r 


2(C, — Cy) — ¢; 


Notes. *\ndependent and identical multivariate normal distribution. **Independent. ***Arbitrary multivariate normal distri- 


bution. ****Arbitrary multivariate distribution. 


it must be based on an economic rationale: The ware- 
house is open to lower the expected cost from leftover 
inventory, but at the expense of incurring incremental 
transshipment costs. Hence, this tradeoff is justified 
when the penalty of lost sales is high enough. Specif- 
ically, using the K-K-T conditions from Table 1 it can 
be shown that it is optimal for the warehouse to be 
used when 


0 < cr/(Cr +¢), — Cy) 
= F,(S;) < F(S*) 
= (GC) (eo Cr Cp) =, (6, ,,) 
aa F(S") 


The significance of (3.1) is that not only does the 
retailer stock less because an economic backup is 
available, but total inventory in the system is lower 
than would be the case if the retailer acted indepen- 
dently by stocking SX’. It is as if the imperative to 
reduce costs is paramount, even though it is at the 
expense of serving fewer customers. 

Not only does the case n = 1 show that there can be 
a pure economic rationale, it also leads to the conclu- 
sion that the simple threshold policy is also optimal 
when there are n retailers whose demand is perfectly 
correlated, because in that case rescaling reduces the 
problem to a one retailer model. 


(3.1) 


COROLLARY 4. If demand at each retailer is perfectly 
positively correlated, then 


WO 
So 


>0@c,>Cy and Cp>C,Cr/(c, —Cy). 


As in the discussion after Theorem 5, Corollary 4 
suggests that the effect on the optimal policy of 
increasing n will be in the opposite direction to the 
effect of increasing variability. As we will discuss later 
this is indeed the case. 


In the remainder of this section we will consider the 
case n > 1. Because it is crucial to the analysis, we will 
focus on characterizing the functional form of © as 
summarized in Table 3. Complementing Corollary 1 
and Corollary 2, these properties lead to establish- 
ing the structure of the optimal policy in terms of 
the parameter c,, and, the sensitivity of this policy 
to the variability in the demand process and system 
size in terms of n, the number of retailers. Finally 
we will conclude with the implications of our results 
for describing the optimal inventory polices in these 
systems. 

To proceed, observe that by Theorem 1.1, 0>0 
when c, < cr. Moreover, using the K-K-T condi- 
tions from Table 1, it can be shown that lim, _,,, 0 
approaches —(c, — cy)/c, and —(c, — Cy)/(Cy — ¢, + 
c, —Cr) for the RF and WE systems respectively and 0 
for the WO system. Additional analysis reveals that 
this limit is reached from below for the WO system 
if and only if —(c, — c/n) is negative. This yields the 
necessary and sufficient conditions in the second row 
of Table 3 that assure that the limit is nonpositive so 
that the warehouse is open for large enough c,. The 
final step is to show that, as for n > 1, © is quasi- 
convex in c,. Conditions under which quasiconvex- 
ity is established are summarized in the third row of 
Table 3. These conditions are quite robust and include 
the desirable case of arbitrarily distributed normal 
demand. These three observations lead to establish- 
ing the structure of the optimal policy in terms of 
the parameter cp. To proceed, let 1(n,2%,) denote a 
nonnegative number, if it exists, that satisfies the 
following: 


So > 0 o Cp> i(n, a) 


i.e., 1 describes a single threshold policy in terms 
of Cp, if one exists, for the warehouse to carry positive 
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stock. Then 


THEOREM 9. Suppose © is qausiconvex in c,. Then, 
(i) For the WO system, 1(n, %,) < 00 C, > C/N. 
(ii) For the RF system, 1(n, %,,) < 00 C, > Cy. 
(iii) For the WF system, 
(a) if Cc, > Cy, then i(n, X,) < ~; 
(b) if c, < Cy, then, either 1(n, &,,) = 00 or there exist 
nonnegative numbers 1,(n,%,) and 1,(n,%,) such that 


S > 0S0< <0o,<b<om. 


Part (i) of Theorem 9 establishes that if the ware- 
house is open as c, approaches oo then a finite single 
threshold policy is optimal. Otherwise the threshold 
is set at oo for the WO and RF system so that S, = 0; 
however, the quasiconvexity of Oy, and a positive 
threshold limit yield a two-threshold policy. This is 
confirmed by numerical examples that show that the 
warehouse is closed for small and large values of c, 
but open for intermediate values, while other numer- 
ical examples show that the warehouse is always 
closed. 

Using the relationship between when a warehouse 
is used and C, established in Theorem 9, we can refine 
Corollaries 1 and 2 to yield the following: 

CoROLLARY 5. When © is qausiconvex in c 
increases: 

(i) If (C3) holds and c, < cy/n < Cy, the optimal con- 
figuration evolves as either NP/RO or NP/RO/WF/RO. 

(ii) If (C3) holds and cy,/n < c, < cy , the optimal con- 
figuration evolves as either NP/RO or NP/RO/WF/RO 
or NP/WO/WE/RO. 

(iii) If (C3) holds and c,, > cy, the optimal configuration 
evolves as either NP/RO/WE or NP/WO/WE. 

(iv) If (C3) is violated and c, < cy , the optimal config- 
uration evolves as NP/RO. 

(v) If (C3) is violated and c, > cy, the optimal configu- 
ration evolves as NP/RO/RE. 


Part (iv) identifies the unique case when backup 
response has only one level. Parts (iii) and (v) con- 
sider the cases when the response is hierarchical in 
the sense that initially the retailer operates in isola- 
tion, then adds the first level of backup response fol- 
lowed by the second level. Moreover, in these cases, 
according to Theorem 9 these policies have the intu- 
itive appeal that once a warehouse is open for values 
of c, above a critical threshold, it remains open for 
all greater values. However, when c; < cy parts (ii) 
and (iv) show that this intuition has to be qualified 
since there can be instances when the dual response 
WF policy reverts to a single level response. 

Now that we have fully established, under the 
conditions of Theorem 9, how the optimal policy 
varies with increasing Cp, We examine the sensitiv- 
ity of the optimal policy to changes in variability 
and size. Specifically, we examine how 1 behaves 
when the warehouse is open. Applying Theorem 6 as 


op AS C, 
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summarized in row 4 of Table 3, we can establish that 
as variability increases © is nondecreasing for the WO 
system, and it is nondecreasing for the RF and WF 
systems when c, <c,, In these cases, if the warehouse 
is closed for a given set of economic parameters, it 
remains closed as variability increases. However, it is 
nonincreasing for the RF system when c, > c, so that 
once the warehouse is open it remains open as vari- 
ability increases. In light of Corollary 4, it thus follows 
that as variability increases, 1 increases to 1(1) > c, for 
the WO system. In contrast, when c, > c,, as vari- 
ability increases 1 decreases to 1(1) for the RF system. 
Hence, a continuity argument yields the following 
powerful result. 


THEOREM 10. In the RF system the warehouse is closed 
when c, < 1(1) = CyCr/(Cy — Cy)- 


While Theorem 10 consistently with Theorem 8 
establishes that in the RF system, 5S) = 0 if OS 1(1), 
it does not reveal how 1 varies with n. Hence, apply- 
ing Theorem 6 as summarized in row 5 of Table 3, 
when the demand distribution is of the class IN, 
Or is monotonically increasing in n when c, > cy, sO 
indeed 1 is increasing in n. And as expected, Oyo is 
indeed decreasing in n, so that 1 is decreasing in n. 

As summarized in the last row of Table 3, a direct 
argument yields that if the demand process is inde- 
pendent at each retailer, 1 approaches cr so that even- 
tually the warehouse is open for all WO systems for 
large enough n. While we are unable to establish such 
a result for the RF and WF systems, the asymptot- 
ical results of Theorem 7 provide an approach for 
establishing such results for the RF and WF sys- 
tems. As summarized in the last row of Table 3, as 
n approaches oo, for the RF system there is a lim- 
ited range of parameters for which the warehouse will 
be open if c, is above a critical threshold which is 
larger than 7(1). However, for the WF system, as in the 
WO system, the warehouse would be open for large 
enough n for all c, > cr if either (a) the transportation 
costs are such that the triangular inequality (c; > c,/2) 
is violated or (b) the triangular inequality is satisfied 
and the parameters fall in a limited range. For another 
perspective on these results, consider their implica- 
tions for Corollaries 1 and 5. Essentially, except for 
a limited range of parameters, the implication is that 
when the choice narrows to NP/RO and NP/RO/RE, 
for sufficiently large 1 NP/RO is preferred. Simi- 
larly, NP/RO is preferred when the choice reduces to 
NP/RO/WE and NP/RO. However, when WO is in 
the mix as in Corollary 1(b) and 5(iii), NP/WO/WE is 
preferred. In other words, in large systems with pool- 
ing, we should tend to see RO systems except under 
conditions of Corollary 1(b), in which case we should 
tend to see the WO system for moderate values of c, 
and the WF system for higher values of c,,. 
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Table 4 Characterization of Inventory 
6) RF WF 
S 2 Semen Yes (A) Yes (A) Yes (identical distribution) 


SG rina No Yes (A) No 


Having established that under robust conditions 
on the demand distribution there is always a set of 
parameters under which it is optimal to carry stock 
at the warehouse, we present the final set of ana- 
lytical results that establish how the optimal use of 
warehouse-based stock influences the stocking deci- 
sions at the retail level. For n > 1, as summarized in 
Table 4, using the K-K-T conditions yields that in the 
RF system, each retailer stocks less than in the bench- 
mark RO system but total inventory is higher. Hence, 
in contrast to the case n = 1, a larger proportion of 
customer demand is met. Similarly, in the WO system 
each retailer stocks less than in its benchmark NP sys- 
tem while numerical examples show that total system 
stock can be higher or lower. In contrast, in all numer- 
ical examples for the WF system that we computed, 
the total system stock was always higher than in the 
benchmark RO systems. And in all problem instances 
each retailer stocked less than in the benchmark RO 
system, a result that we have been able to establish 
analytically only for the case when all retailers are 
identical. 

Using the technical results of §1 on the structure 
of the optimal policy and the extensive results of §2 
for the RO system, we have presented a compre- 
hensive analysis of WO, RF, and WF systems. We 
have established the form of the optimal policy and 
how it is influenced by variability and system size. 
We have also traced the impact of optimally using a 
warehouse-based stock on inventory allocation at the 
retail and the system levels. 

Broadly speaking, these results derived from 
Tables 3 and 4 are very robust and quite complete for 
the WO and RF system. While the results for the WF 
system are at the same level of robustness in terms of 
demand distributions, they are not as extensive. We 
attribute the contrast to two attributes of the neces- 
sary and sufficient conditions presented in Table 2. 
The optimality of the equal fractile rule for each of 
the WO and RF systems obtained from Table 1 is 
used to characterize either the joint distribution or the 
marginal distribution of demand. In contrast, for the 
WE system, the interaction of the joint and marginal 
distributions of demand must be untangled but with- 
out the benefit of the equal fractile rule. 


4. Concluding Remarks 
Motivated by recent reports of improved customer 
service facilitated by the use of transshipment to pool 
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inventory, we developed a formal model that focused 
on the role of transshipment in a multiple retailer net- 
work. In §1, we determined that for the dynamics of 
our model, the optimal policy precludes the use of 
partial fulfillment of backordered demand. As much 
as possible of backordered demand is filled or none 
at all. Moreover, the transshipment protocol when it 
is optimal is represented by exactly one of four sys- 
tems. The easily computed conditions that determine 
the optimal policy are summarized in Theorem 2. The 
final test requires a simple computation which at most 
requires the solution of a functional equation for a 
unique value. Moreover, Corollary 1 and the unified 
analysis of the three warehouse-based systems shed 
additional light on the form of the optimal policy. 

If the cost of holding inventory is higher at the retail 
level than at the warehouse, then our results suggest 
that for large enough penalty cost of stocking out it 
will be economical to use some sort of pooling pol- 
icy. Most significantly, the WO system is only opti- 
mal for relatively intermediate values of the penalty 
cost under a restricted set of parameters. Under such 
parameters, from a practical perspective, it is easy to 
justify the use of the warehouse in all but small sys- 
tems. This is so because the threshold for c, is decreas- 
ing in the number of retailers. In contrast, when the 
choice is between the RF and RO system, the balance 
will tend to tilt away from the RF system as the num- 
ber of retailers increases. Intuitively, one can explain 
this result by noting that this is so since the ware- 
house is the provider of last resort among increasing 
number of backup providers. Finally, the case when 
the choice is between WF and RO systems is more 
ambiguous. If the transportation costs are such that 
the triangular inequality is violated, improbable in 
most practical settings, then the economic benefit of 
shipping from the warehouse will propel the use of 
the warehouse, which is now the first of many back- 
ups. However, if the triangular inequality is satisfied 
there is only a limited range of parameters under 
which the warehouse would be used. In particular, 
the WF system would be optimal for large enough n 
(since the WO system would be optimal for interme- 
diate values of c, such that transshipment at the retail 
level is not economical). Hence, in summary, for a 
small system with intermediate values of penalty cost 
the value of pooling from a warehouse is limited. 

The single-period models considered here are nat- 
urally related to periodic review settings, a multi- 
period setting for modeling inventory control systems 
for durable and consumer goods. For example, our 
model applies directly to multiperiod settings by gen- 
eralizing the model of Robinson (1990) if like him 
we assume that the lead time for a restock order is 
zero. This multiperiod model can be embedded into 
infinite-horizon variants as the building block for the 
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approximate single-cycle analysis like that in Eppen 
and Schrage (1981) and Schwarz (1989) for variants 
of the WO system. We refer the reader to Chapter 5 
of Wee (2000) for the extension of the single-period 
model considered here to such multiperiod models. 

The results of our paper would predict that in 
these periodic models the optimal transshipment pol- 
icy would exclude partial transshipments; however, 
these results have to be qualified when more com- 
plex dynamics are involved. For example, Anupindi 
et al. (2001) consider a system in which retailers 
are independent sellers of a manufacturer’s prod- 
uct; in such models because of strategic interactions 
between different parties the manufacturer may facil- 
itate storing some inventory at a central warehouse 
to fill backordered demand. Moreover, the indus- 
trial settings of after-sales service networks like those 
described by Narus and Anderson (1996) use the 
warehouse for such shipments. However, the dynam- 
ics of those systems are fundamentally different: 
Inventory is reviewed using continuous time models 
and the warehouse plays the dual role of filling reg- 
ular replenishments and emergency shipments. Some 
analytical insight into this difficult model is provided 
by the heuristic approach presented in Grahovac and 
Chakravarty (2001). Finally, another perspective on 
the complexity of these continuous review models can 
be gleaned from a series of papers by Axsater (2003). 
In a recent paper in this series (and related references 
therein), a continuous time version of the RO system 
is considered. Because replenishments have a positive 
lead time, in the same spirit as Archibald et al. (1997), 
a transshipment policy that may be interpreted as fea- 
turing partial pooling is proposed and shown to be 
effective. 
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Appendix 

DERIVATION OF NSywo, NSpgr, AND NSwr. We shall only 
consider NSwo since the other two cases are similar. Notice 
that Sj’° =0 if and only if there exist corresponding SW°, 
i=1,...,n and A;,i=0,...,n such that they satisfy the 
K-K-T conditions. Let 55° =0. Then, observe that H,,(0) = 
P(N,{D; < S¥°}) = F(S#O,..., 5°) and dE(Lw)/AS;l5-0 = 
P({D; > S%°}.M {Dy < 0}) =0. Thus, from the K-K-T condi- 
tions, Ay = (cy +c, — Cr)H,,(0) + cr — c,, E (SY) = c,/ (Cy + 
c,) = F(SN”) and A; =0, i=1,...,n. 

Clearly, the above solution satisfies (1.11)-(1.14). Hence, 
the K-K-T conditions are satisfied if and only if Ay > 0. Note 
that from above, Ag = (Cy +c, — Cr) F(SS”,..., SN?) +7 — Cy. 
Hence, Ap > 0 if and only if (cy +c, —cr)F(SY",..., SNP) + 
Cr — C, = 0. Rearranging the last inequality, we have (NSwo) 
and the result follows. 

PROOF OF LEMMA 2. Because the function E(Ly) is 
integrable and has a bounded derivative, it satisfies the 


Management Science 51(10), pp. 1519-1533, © 2005 INFORMS 


Lipschits condition of order one and hence the expec- 
tation and derivative in dE(Lw)/dS, can be interchanged 
(see Glasserman 1994). Hence, dE(Ly)/0S, = E(dLy/0S,) = 
E(dLw/dDp x 0Dp/0S,). Let 1,,, be the indicator function of 
the event w, ie., 1,,, =1 if @ is true and 1,,, =0 otherwise. 
Using the indicator function, we have 


0E(Lw)/0S, = E((—1)1ip,55, x (~1)1)p,<5))) 
= P({D, > 51} {Do < So}). 


PROOF OF COROLLARY 3. We will focus on (iii) since (ii) 
can be proven in a similar way while (i) is a direct con- 
sequence of (ii) and (iii). Suppose C,>c, but on the con- 
trary, F,(SK°) < 1/2. This implies S®° < y; Vi and hence, 
G,(S®°) < 1/2 < ¢,/(c, + ¢,), a contradiction to Theo- 
rem 4(iii). Hence, F,(S§°) > 1/2 and S®° > w; Vi. Let z, = 
(SFO — ,)/o}. 

Let ® denote the c.d.f. of the standard normal distribu- 
tion. Then, 


6. (SRO) — o(D18° -u)/ (x o} +Eo) ) 
; 1 t,] 


-0(sE0/(Ee+En)") 


> @(z,) (because z, > 0 and the Schwarz 
inequality holds) 
= A (5°). 


Hence, from Theorem 4, we conclude that G,(S%°) > 
F, (SN?) > F,(SR°) > 1/2. 

PROOF OF THEOREM 7. We will focus on (ii) since (i) can be 
proven similarly. Suppose c, > c;, and let ® denote the c.d.f. 
of the standard normal distribution. Because demand at 
each retailer is identically distributed, each retailer has the 
same optimal stocking level so that G,,(S®°) = ®(n(SR° — 
Mi)//noz) = ®(z,./n). Hence, as n — oo, if z, does not 
decrease to 0, then G,(S%°) increases to 1 and from (2.1), 
F,(SK°) decreases to (c, — c,)/c;. This is possible only if 
(c,—c,)/c, => 1/2 since F,(SK°) > 1/2. Otherwise, z, decreases 
to 0, ie., F.(S®°) decreases to 1/2 and G,(S®°) increases to 
(c, — €,/2)/(Cy +c, — ¢). 

PROOF OF THEOREM 10. 

Case 1 c, <c,. From Corollary 3(i) and (ii), R(S%°) < 1/2. 
Hence, (NSpr) is satisfied if the right-hand side of (A3) 
is greater than 1/2, which reduces to c,(2cy +c, — 2c,) + 
2C,Cr — CyCp — CyCy = 0. 

Subcase 1 2c, +c, — 2c, <0. Then, 

Cy(2cy + Cy — 2Cy) + 2C, Cp — CyCp — CyCy 
> C,(2cy + C, — 2c,) + 2cy,C7 — Cp — CC Hy 
= (2cp, — C;)(Cy — Cy + Cr) 
5 2 Oy G, => 2c; = C). 
Subcase 2 2cy+c,—2c, > 0. Then, 


>0O (since c,>Cc 


Cy(2Cuy + Cy — 2C,) + 2€, Cp — CyCp — CyCy 
> (Cp — Cy) (2€y + C, — 2c,) + 2€,C7 — CYC — CCH 
= 21 (Cy — Cy + Cr — C;) 


> 0. 
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Case 2 Cy < Cy < CyCr/(C, — Cy). Since cy, < Cy, from Corol- 
lary 3(iii), F(S°) < R(SS*) = c,/(c, + c,) so that (A3) 
reduces (¢, +€), —C;)(CyCH — Cy, +CrC},) = 0. Because c, + Cy, — 
c, = 0, (A3) holds if c,cy —¢,c, + cre, = 0. If c, < cy then 
the last inequality clearly holds, and if c, > cy then the last 
inequality reduces to c, < CyCr/(Cy — Cy). 

DERIVATION OF TABLE 3, Row 2. For the WO system, note 
that lig, 00 @Owo =1-—1=0. Moreover, 


_ O@Mwo  ,. ( r OF 1 ) Ch CH 
lim = lim — « —— _| ——_ — ————... 
Cpe OC, = pre 2. 08;  fAS)) (Gp+en)? (cp, +¢ey—-er)? 
Note that 
lim ae x : 
1 a 
ce AS; — fi(S)) 
= = Fis, sep Up awe, ay de 
Hence, 
lim aed = lim ead eee ~ ————#__. 
cpr OC, Cpe (Cr — Cy) cocoa Cr) 
Thus, lim, ,,, 9Owo/dc, = 0* if and only if nc, > cy and the 


result follows. Similarly, for the RF and WF systems, the 
theorem follows by noting that 


C;—C Cy —C 
lim O.4e1—— ee 
is Ce = GC, = ty Cy =O, FG — Ce 
DERIVATION OF TABLE 3, Row 3. 
WO System: @wo 1s quasiconvex in C,. 
We shall show that 0° @y-/0c5|s0 /acy=0 = 9 


Iwo _ SOF 7 0S; CH 
aC, 09; OC. (Cy Cae Cp) 
_ OF 1 Ch _ Cy 
i=l 0S; FAS) )* (c, 7 C;,)? ic + Cy cr)? 
C 0S 1 C 
F(S, =—- = — = —— x —_ > 0) 
(HO= ee Gea 
fh Cy 
7 (c, + C),)? (c, sy CH — Cr)? 
where A = {>-j_, OF /0S; x 1/f;(S,)}. 
la) 
WO _ “ 9A Ch ’ 2Ch x 2Cy . 
dc; (Coe C,)2 (Coc)? Ac, ey — Cy)” 
POwo 
dc; 4/Acy=0 
JA Ch, 


x ec 
O07. Key Oy)? 


+2] ot 7 ee 
(Cp+¢y—Ccr)?  (C,+¢y — Cr)*(cy + Cy) 


0®) Cc C 
(since NO <0 Ax th =H 
ac, (CG) ACeey— cy) 
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0A C 2c 
0G, (CG) (Co Ca = op) 
Peery 
Cytcy—Cp Cy tty | 
Since 1/(c, + Cy — Cr) — 1/(Cy + ¢)) = (Ch + Cr — CH)/(Cp + Cn) - 


(c, +c — Cr) = 0, it remains to show that dA/dc, > 0. Note 
that A = >°7_, A; where 


A= - [op f- fi lO)... px 


Xn 1X; = 5;) I] Ax, 
k#i 


Xj- —1/ a 


and f,""' is the conditional distribution of the n —1 random 
variables x,,...,Xj-1,Xj41,+-+,%X, given that x; = S;. We now 
show that 0A,/dc, > 0. We first partition 2, as follows: 


— 07 7 
Xa Ly 
Then, 
So Sy 
A= | _ fA Gar, 


where f”"' has a multivariate normal distribution with 
mean yp and covariance matrix K: 


X, |X, =5S,)dx,---dx 


ns 


, PinFn) (Sy ~~ by)/O1; and 


(see, for example, Graybill 1976). 


b= (My, eg bly) (Pp Oop +s: 
K =X» — Lot 212/07 


Applying a change of variable y; = x; — uw; — p4;0;(S 
and the identity (5; — w;)/o; = (S; 


1 Hy)/O4 
— ,)/G, for all 1, we have 


0 (1—p12)(Sy —1)/O%74 On (1—Pin)(S) —H1)/0%1 
A= | of Di (Yor +61 Yn) 4Yp*° Yn, 
(A2) 


where ¢,; is the standard multivariate normal distribution 
with covariance matrix K. From (A2), we see that 0A,/dc, = 
dA,/0S, x 0S,/dc, = 0. Similarly, dA;/dc, > 0 fori =2,...,n 
so that dA/dc, > 0. 

RF System: Opp is decreasing in c,. 

The result follows since 


OOrrp — —Cy(Ch+Cr—Cy—C;) OF, 


dc, (Cy + — Cr)? dc 


Differentiating (1.8) with respect to c,, we have 


aS; _ aS fils) 


- | (A5) 
dc, dc, — f(S)) 
We shall show that @°@ywr/0C5| 56 /acy=0 = 9. 

0®) “OF 9S, C 0S 
ae ee -(—+ ) x fi(S1) x 5 
p eee Oy CH Cyt Cy — Cr Cy 
"OF 1 ¢, x aS, 

= |e x pe I Al x 
» 0S; FAS) (Cy —C,+C;—Cr) *f(G) dC, 


Hence, F Owe /9Chla@/ac,-=0 = dA/dc,, where A = )°7_, OF / 
0S; x 1/f;,(S;). Since the equal fractile rule continues to hold, 
ie., (S; — w,;)/o; = (S; — p,)/o, for all 1, Owele,=cr > 0 and 
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Owele,=c <0, the rest of the proof follows that for the WO 
system. 

DERIVATION OF TABLE 3, Row 4. The result for the WO 
system follows from (NSwo) and the observations that SN? 
is invariant in %, while F(S,,...,5,) is lexicographically 
increasing in &, for fixed S;,’s (see Plackett 1954), while the 
result for the RF system follows from a direct application 
of Theorem 5. Hence, it remains to show that @w, is lexico- 
graphically increasing in X, when cp <c, for any arbitrary 
multivariate normal distribution. Let 


crf (Sf°) = er 


0, (SE, ... 
Cp i, Pp ee 


yon? = (Soc .7 SEO) = 
when the covariance matrix is 2. Suppose X, and &, are 
such that o;, > Ti for all i, j; it suffices to show that ©, > @,. 
Let 2 4(,) =AX_ + (1 —A)z,. Then, 


JO,)/9A 
= 0F/0A + _ dF /dS° x aSKO/AA 


—fi (ST es) (cx —C,+¢,—Cr) x aSFO/dA 


ay ee Ge ~ Cy +¢,— Cr) — )\(AF /aSF°)/f, (s*°)| 


x fi(SF°) x [as®°aal 


(since 0S%°/dA <0 when C= C,). 


Claim. For c, <¢), f;(S$°) > n. 


For convenience, let f* denote f(x,,.. 
Note that 


pl agate as 


OF ph phn 
ag 1 LAT 
coef fT] ax 


7 k#i 

es a ees es 

e ae we me oe 
Lo). Pant {fof {Pla 


n—1 
+f ff pies 


k#i 


flS) = f 


Ee EE Le te) 


k=0,kAn-1 


+f [pT 


7 kAi 


IV 


a 5 
n i, vee i fT] ax, 
oe = k#i 
(since S, < @, for all k when c, < c,). 
Hence, >7;(0F /dS®°)/f,(S®°) < 1. Since F/dA = 0 (see Plack- 


ett 1954) and c,/(cy —c, +c; —Cr) = 1, 00,/0A = 0. In partic- 
ular, 0, = Oo(1) = 0.00) = @.. 
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DERIVATION OF TABLE 3, Row 5. The result for the 
WO system follows from (NSwo) and the observation that 
F(SN?,..., SN?) is nonincreasing in n, while that for the 
RF system follows from a direct application of Theorem 6. 
Finally, the proof of the result for the WF system is similar 
to that for row 3 for the WF system. 

DERIVATION OF TABLE 3, Row 6. 

WO system 

The result follows by observing that when demand is inde- 
pendent, F(SN?’,..., 5%?) = (F,(SN’))" which is decreasing 
in n and converges to 0 as n diverges to oo. 

RF system 

Since S¥ =0 for ¢, < c,Cr/(c, — Cy), we need only to con- 
sider the case when c, > c,Cr/(C, — Cy) = C,. From Theo- 
rem 7, we have 


P 


Cyl Cyr, =C 
CNG) if c,>2c,, 
€;(c, +€y—Cr) 

lim Ogr = 

are C207 FO — 2G) = Cr (C20, ) Cae 


if c, >2c,. 
2¢;(Cy + Cy — Cr) ‘ 


Note that the first limit is always positive. The numerator 
in the second limit is nonpositive if and only if 


2CyCr — C,(Cp + Cy) 


P= “2c, —e,) —G, and 2(c, —Cy)—c, > 0. 


(Note that if 2(c, — cy) — cr <0, then c, < [2c,cr — ¢,(cr + 
Cr) /[2(c, — CH) — Cr] < Cyer/(Ch — Cy)-) The theorem then 
follows by rearranging the second condition and apply- 


ing (PC2). 
WE system 
For convenience, let B = cy —c;, +c;—Ccy. From Theorem 7, 
(cr — c,/2)/B if cy, > Cy > ¢,/2; 
Cr—C <Q) ifc,>c,andc,/2>Cc,; 
fim Og. = {OT VB (<0). if eye and c4/2> ¢, 
ae (cr — ¢,/2)/B if Cy > Cy > C,/2; 


(Cr —c, + ¢),)/B 


The result then follows by applying the appropriate limits to 

the cases cy <c,/2 (with subcases c, < Cr, Cr < Cc, <c,/2 and 

c,/2 <c,) and cr > c,/2 (with subcases, c; <c,/2, c,;/2 <c, < 

Cr, and cr <c,,). Note that for the subcase c, < c,/2 < cy, the 

warehouse is open if and only if c, > c),, +c;. To complete the 

proof, note that c, > c, +c; and c,/2 < cr together imply that 

C, < ¢,/2 (<cy) so that this last condition is automatically 

satisfied and hence is not included in the last cell of Table 3. 
PROOF OF TABLE 4, Row 2. 

WO System 

Suppose on the contrary S/‘° > SN? for some k, then 

F(S°°) > F(SY") = c,/(c) + cy). The last inequality 

and (1.12) imply that dE(Ly,)/0S; < 0, a contradiction. 

RF System 

Observe that because of the equal fractile rule, if 4k such 

that SKF > SRO, then SKF > SRO > 0 Vi. The proof is then 

complete by comparing (1.8) with (1.22). 

WF System 

Comparing (1.8) with (1.17), yields that there exists an 1 

such that S/F < S®°. (Otherwise, the left-hand side of (1.8) 

is positive.) The result then follows since all retailers are 

identical. 


Hf Cy > Cp aNd Cy/ 2 > Cy. 
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DERIVATION OF TABLE 4, Row 3. We only need to prove 
the result for the RF system. Suppose on the contrary that 
SRF < SRO, Then, G,,(7_, SR) < G,(SR*) < G,,(S®°). Substi- 
tuting the above two inequalities into (1.22), we have 


(c, +6, —¢,)G,(S*°) +.¢,5(SF°) —cp <0, i=1,...,n, 


thereby contradicting (1.8). 
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